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We present two experiments that provide information on spin anisotropy and the magnetic struc- 
ture of Na^CoC^. First, we report low-energy neutron inelastic scattering measurements of the 
zone-center magnetic excitations in the magnetically ordered phase of Nao.7sCo02. The energy 
spectra suggest the existence of two gaps, and are very well fitted by a spin-wave model with both 
in-plane and out-of-plane anisotropy terms. The gap energies decrease with increasing tempera- 
ture and both gaps are found to have closed when the temperature exceeds the magnetic ordering 
temperature T m » 22 K. Secondly, we present neutron diffraction studies of Nao.8sCo02 with a 
magnetic field applied approximately parallel to the c axis. For fields in excess of ~ 8 T a mag- 
netic Bragg peak was observed at the (0, 0, 3) position in reciprocal space. We interpret this as a 
spin-flop transition of the A-type antiferromagnetic structure, and we show that the spin-flop field 
is consistent with the size of the anisotropy gap. 

PACS numbers: 75.40.Gb, 74.20.Mn, 75.30.Kz, 75.25.+Z 
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I. INTRODUCTION 

Since the discovery of superconductivity in water- 
intercalated Na^CoO^i much of the discussion has fo- 
cussed on the mechanism of superconductivity. There 
is strong support from experiment^ and theory^ for an 
unconventional pairing state, the origin of which derives 
from the triangular lattice of the Co ions and the exis- 
tence of strong spin and charge fluctuations. However, 
details of the pairing state are far from resolved^ The 
phase diagram of Na x Co02 shows two magnetically or- 
dered phases, one at x s=s 0.5 and the other in the range 
x » 0.7 — 0.95. The presence of these suggests that mag- 
netic correlations may play an important role in the for- 
mation of superconductivity in hydrated Na x Co02, as is 
believed to be the case for the superconducting cuprates. 
This possibility provides a strong incentive for character- 
izing the magnetic order and excitations of Na^CoO^. 

This paper is concerned with the weakly magnetic 
phase found for x w 0.7 — 0.95 at temperatures be- 
low T m w 22 Km& Recent neutron scattering studies 
of Nao^CoOa^ and Nao^CoO^ have established that 
the magnetic order and dynamics are consistent with an 
A-type antiferromagnetic structure, shown in Fig.^ and 
that the magnetic interactions are three-dimensional de- 
spite the two-dimensional character of the crystal lat- 
tice and electronic structure. Reasons for the relatively 
strong c-axis magnetic coupling have been discussed re- 
cently by Johannes et alA2> Surprisingly, although the 
magnetic correlations are rather strong the ordered mo- 
ment is only 0.1-0.2 fj, B M^ 

The first aim of this work was to extend the measure- 
ments of the spin-wave dispersion in Nao.7sCo02 down to 
lower energies where we might gain important informa- 
tion on the magnetic ground state, such as whether itiner- 
ant effects are important. Our previous measurements, 8 
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FIG. 1: (color online) (a) Crystal structure of Na;cCo02. (b) 
The A-type antiferromagnetic structure on which the spin- 
wave model is based, showing the two exchange constants J a t> 
and J c . (c) The a-b planes, showing the orientation of the 
oxygen tetrahedron around a central Co ion. 



as well as those by Bayrakci et al. on Nao.82Co02r found 
evidence of a small excitation gap at the antiferromag- 
netic zone center. Here we present a detailed study of 
the spectral lineshape in this low-energy region which 
not only confirms that the magnetic excitation spectrum 
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is gapped, but also provides evidence for the existence of 
two gaps. Both these gaps are found to close at the bulk 
magnetic transition temperature T m «22 K. 

The second part of the paper presents a novel way 
to investigate the magnetic structure of Na x Co02 (x w 
0.7 — 0.95) and also to gain information about spin 
anisotropy. In the accepted A-type antiferromagnetic 
structure (Fig. ^>) the moments are ferromagnetically 
aligned within the layers and stacked antiferromagneti- 
cally along the c axisiS^ The ordering wavevector for this 
structure is (0,0, 1) Strong spin-wave scattering is ob- 
served emerging from (0, 0, 1) positions with odd I, which 
are zone centers for the A-type antiferromagnetic order, 
but no magnetic Bragg peaks have been observed at these 
positions using neutrons.- On this basis it was deduced 
that the ordered moments point along the c direction 
since neutrons scatter from the component of the mo- 
ments perpendicular to the scattering vector. This mo- 
ment direction is consistent with what has been inferred 
from the uniform susceptibility and from //SR data. 12 
Moreover, Bayrakci et al. did succeed recently in ob- 
serving magnetic Bragg reflections at a few (h, k, I) posi- 
tions with h, k ^ and odd I using polarized neutrons^ 
Again, these are consistent with the accepted magnetic 
structure. Polarized neutrons were required because the 
ordered moment is small and strong non-magnetic scat- 
tering is observed at all positions where magnetic Bragg 
peaks are expected. Up to now no magnetic Bragg peaks 
have been observed with unpolarized neutrons. 

The experiment we describe here was originally de- 
signed to confirm the proposed magnetic structure by a 
method that employs unpolarized neutrons but avoids 
the problem of having to separate the weak magnetic 
scattering from the strong non-magnetic background sig- 
nal. Our approach was motivated by measurements of 
the magnetization of Nao.85Co02 in applied fields up to 
14 T by Luo et al.m& The magnetization data with H [| c 
show a clear anomaly at 8 T at low temperatures, with no 
such transition seen for H _L c. The authors interpreted 
this as a spin-flop transition in which the ordered mo- 
ments rotate by ~90 degrees while preserving the A-type 
antiferromagnetic arrangement. After the transition the 
spins lie approximately in the hexagonal plane, but the 
magnetic structure has a small ferromagnetic component 
along the c axis. Assuming this explanation to be correct 
we induced the spin-flop transition in a neutron diffrac- 
tion experiment and searched for magnetic Bragg peaks 
along Q = (0, 0, 1), since now the ordered moment should 
be perpendicular to the scattering wavevector and should 
scatter neutrons. Our experimental results are in excel- 
lent accord with the predicted behavior. Encouraged by 
this we go on to show that the size of the spin gap in the 
magnetic excitation spectrum is in agreement with the 
observed spin-flop field. This quantitative analysis pro- 
vides the link between the static and dynamic magnetic 
properties explored in this paper. 

The remainder of the paper is organized as follows. In- 
elastic neutron measurements probing the low-energy re- 



gion of the spin- wave dispersion are presented in the next 
section, along with the analysis of these results. Diffrac- 
tion measurements made with an external magnetic field 
applied parallel to the c axis are presented and analyzed 
in Sec. III. Section IV contains a discussion of the results, 
and the conclusions are presented in Sec. V. 

II. INELASTIC MEASUREMENTS 
A. Experimental Details 

Inelastic neutron measurements were performed on the 
cold-neutron triple-axis spectrometer IN14 at the Insti- 
tut Laue-Langevin. This instrument was chosen to allow 
investigation of the excitations in Nao.7sCo02 at lower 
energies than previously studied^ We employed a py- 
rolytic graphite (PG) (002) monochromator and a PG 
(002) analyzer, which were curved vertically and hori- 
zontally respectively to maximize the count rate. The 
majority of measurements were made with a fixed final 
energy of Ef = 4 meV. A Beryllium filter was placed in 
the scattered beam to suppress higher-order harmonics. 

The inelastic neutron measurements were performed 
on the same crystal of Na .7 5 CoO2 as used for our previ- 
ous experiments which investigated the spin fluctuations 
in this compound at higher energies^ The single crystal 
of mass ~ 1.5 g was mounted on a copper mount and 
aligned to allow measurements to be made within the 
(100)-(001) scattering plane. Details of the growth and 
mounting are given elsewhere&i^ 

Previous measurements revealed strong spin-wave 
scattering around (0,0,1) and (0,0,3), consistent with the 
proposed A-type antiferromagnetic structure, and the in- 
elastic data presented here concentrates on spin waves 
dispersing from the magnetic zone center at (0,0,1), 
where the inelastic scattering is most intense. 

B. Results 

Figures |5] and |31 present examples of inelastic neutron 
scattering data collected on IN14. Each scan was per- 
formed by measuring the intensity of scattered neutrons 
as a function of energy transfer up to ~ 3 meV at the 
wavevector Q = (0,0, 1). Scans were made at ten tem- 
peratures between 1.5 K and 24.4 K. 

Figuresl^ - ^ show the measurements made at low tem- 
perature (T = 1.5 K). The spectrum consists of an in- 
tense peak due to incoherent nuclear elastic scattering 
centered on E = meV, and a broad signal centered 
around 2 meV which is attributed to magnetic scattering 
as the scan cuts through the spin- wave dispersion. There 
is clearly a gap where the intensity falls to background be- 
low ~ 1 meV, revealing that the magnetic excitations in 
Nao.75Co02 are separated from the ordered ground state 
by a clean gap. To determine the non-magnetic scatter- 
ing an energy scan was also made at Q = (0, 0, 1.25). The 
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FIG. 2: (color online) Neutron inelastic scattering from 
Nao.75Co02 measured at T = 1.5 K. (a) Energy scan at 
constant Q = (0,0,1), compared with the same scan at a 
background position Q = (0,0,1.25). (b, c): The data at 
Q = (0,0, 1), fitted (solid line) with (b) a two-mode disper- 
sion (Eqn|2l, and (c) one mode only (Eqn|5]but with E = 0). 
The two peaks under the data in (b) show the contribution 
of each mode to the total intensity, while in (a-c) the dashed 
curve represents the contribution of the incoherent peak and 
background. The shaded peaks show the incoherent peak con- 
tribution scaled down by a factor of 100, as an indication of 
the instrumental resolution, (d) Qz scans with constant en- 
ergy transfers of 0.5 meV, 1.0 meV and 1.6 meV. Data at 
1.0 meV and 1.6 meV have been shifted up by 100 and 200 
counts respectively for clarity. 



scan, which is plotted in Fig. Oa, contains the nuclear in- 
coherent peak together with a small constant background 
signal. Figure [2Ji displays Q; scans performed at three 
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FIG. 3: (color online) (a)-(d) Energy scans at Q = (0, 0, 1) 
at temperatures between T = 1.5 K and 24.4 K. Solid curves 
represent the best fit of Eqn|5]plus an incoherent peak to the 
data. The dashed line shows the contribution of the incoher- 
ent peak and background, while the pink and green curves 
represent the intensities of each mode. 



constant energy transfers of Et = 0.5, 1.0 and 1.6 meV. 
The peak present at higher energies has clearly disap- 
peared at 0.5 meV, confirming that the intensity of the 
spin- wave dispersion really does fall to background in the 
'gap', and the remaining intensity at this point seen in 
Fig.s[5^-c is simply due to the tail of the incoherent peak. 

FigureOlshows the same scan made at three more of the 
temperatures we measured, below and above the mag- 
netic transition temperature T m w 22 K. It appears that 
the magnetic scattering intensity moves lower in energy 
as the temperature increases. Somewhere around 20 K 
the gap seems to disappear, moving into the incoherent 
peak. 
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C. Analysis 



S yy (Q,co)^. where 



To determine whether the gap is in fact decreasing with 
temperature, or whether what we see is simply the mode 
broadening and decreasing in intensity, it is necessary to 
compare the experimental results with a model of the 
excitations. The scan observed at 1.5 K (Fig|2i) is sug- 
gestive of a two-peak lineshape. In order to extend the 
spin-wave model introduced in our previous work 8 to al- 
low two non-degenerate gapped modes at Q = (0, 0, 1) 
we include two anisotropy terms. The Hamiltonian of 
this refined model is then: 

Ti = J a b S ' Sj • Sj' + J c S ' Sj • Sj 

-Dj2(S?f-Ej2[(S?f-(S?) 2 ] , (1) 

i i 

where J a b and J c are intra- and inter-layer exchange 
parameters, respectively, as in our previous work. 8 
The anisotropy constant D quantifies the tendency of 
the spins to lie along the c axis^S while the term 
EJ2i[(Sf) 2 - (Sf) 2 ]i which has two-fold symmetry in 
the plane, is the simplest way to introduce in-plane 
anisotropy. We define x parallel to a, z parallel to c and 
y perpendicular to x and z so as to make a right-handed 
set. From the hexagonal arrangement of the Co ions 
within the a—b layers of this compound one might pro- 
pose a term with hexagonal symmetry for the in-plane 
anisotropy. However, with the quantization direction 
parallel to the c axis a term with hexagonal symmetry 
does not lift the two-fold degeneracy of the spin-wave 
dispersionjSi so for the purposes of this analysis we use 
the Hamiltonian above (Eqn. 

The spin- wave dispersion resulting from this Hamilto- 
nian can be calculated using standard methods, giving 
two modes: 



S xx (Q,lo) = 2S' 



2 {(A Q + D)-(C Q -E)} 



hoj± = 2S^J(A Q + D) 2 -{C Q ±E) 2 , (2) 

where Tiw is the energy transfer, S is the spin (here as- 
sumed to be S — 1/2), Q = (h,k,l) is the wavevector, 
and Aq and Cq are defined in terms of the exchange 
couplings as 

Aq = Jab{cos(2irh) + cos(27rfc) + cos[2tt(/i + k)] - 3} 

+ Jc 

Cq = J c cos(nl). (3) 

The magnitude of the gaps at the magnetic zone center 
are then related to the exchange and anisotropy param- 
eters as follows: 



huj^ = 2Sy/( J c + Df - ( J c ± Ef 



(4) 



Note that if E — only one gap results. 

For the case when Q lies parallel to the ordered mo- 
ment direction, such as at Q = (0, 0, 1) here, the in- 
elastic neutron intensity is proportional to S xx (Q,lo) + 



Q 



xG-(,- Wq )/ 2 (Q)H u )+1] (5) 



S vv (Q,lo) = 2S' 



{(A Q +D)-(C Q +E)} 



Tilo 



Q 



xG+^-c+)/ 2 (Q)[nH + l] , (6) 

where [n(u>) + 1] is the Bose factor, and /(Q) is the form 
factor (which is a constant here since all our energy scans 
were made at one fixed value of Q). G(u> — loq) are 
normalized Gaussian functions which replace the usual 
Delta functions to allow inclusion of intrinsic broadening 
of the two modes. 

The triple-axis spectrometer has a three-dimensional 
ellipsoid-shaped resolution, and therefore does not probe 
the dispersion relation at an infinitely sharp point in re- 
ciprocal space. In order to fit the spin-wave model to 
the experimental data it was therefore necessary to con- 
volute the calculated spectrum (Eqns. OandEJ) with the 
INI 4 spectrometer resolution. This was achieved using 
RESCAL, a set of programs integrated into Matlab which 
calculates the resolution function of the neutron triple- 
axis spectrometer ii& It allows simulation of scans using 
a 4D Monte-Carlo convolution of the resolution function 
with the specified spectrum, and the simulation can then 
be fitted to the data in order to extract the parameters. 

In this way the anisotropy parameters were extracted, 
while fixing the exchange parameters J a b and J c to their 
previous values of —6 meV and 12.2 meV respectively? 8 " 
The relative amplitudes of the two modes were fixed 
by the spin-wave model, with an overall amplitude fit- 
ted, the values for the intrinsic widths of the dispersion 
modes were fitted independently. The incoherent peak 
and background were included as a fixed Voigtian peak 
plus a constant. 

At T = 1.5 K the values for D and E were found to be 
0.096±0.005 and 0.059±0.005 meV, corresponding to two 
modes with gaps of 0.95 ± 0.13 and 1.95 ± 0.06 meV. To 
achieve a good fit the intrinsic widths of the two modes 
were found to be different, 0.37 and 0.74 meV for the 
lower and higher modes. The fitted curve is displayed 
on Fig|2K) with the two lower curves representing the 
contribution of each of the gapped modes. The dashed 
line shows the contribution of the incoherent peak, which 
is also plotted scaled down by a factor of 100 (shaded 
peak) as an indication of the instrumental resolution. 

The model appears to fit the data well. For comparison 
Fig. Et> shows the same data fitted with a 'one-mode' 
dispersion, by fixing the value of E to zero. It is clear 
that the model with two modes fits the data better than 
that with one, yielding a value of x 2 = l-2 compared to 
X 2 =3.7 with E = 0. In Fig. 0] we plot the two dispersion 
modes parallel to Q; (calculated from Eqn. [21 using the 
fitted parameters for D and E), together with the data 
previously measured around Q = (0,0,3). The fitted 
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modes are also in good agreement with the experimental 
data in the Q; direction. 
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FIG. 4: (color online) Filled circles: Magnon dispersion 
parallel to (0,0, 1) centered on (0, 0, x)=(0, 0, 3), measured 
previouslyi- Open square/diamond: Fits to 1.5 K energy scan 
at (0, 0, x) = (0, 0, 1) shown in Fig|5] Solid curves are modes 
calculated from the spin-wave dispersion Eqn|2]witli the con- 
stants given, while dotted lines represent the zone boundaries. 

The fitting procedure was repeated for data at all tem- 
peratures, restricted only by fixing the intrinsic widths of 
the two modes to the values at 1.5 K. Figures |3t~d are 
overplotted with fits to each data set, with the contri- 
butions from each mode as solid lines underneath, and 
the contribution from the incoherent peak denoted by 
the dashed line. The fitted lines provide a reasonable 
description of the data, reproducing the shift of the in- 
tensity towards zero energy with increasing temperature, 
although clearly the lineshapes fit less well as the temper- 
ature increases. The gap energies extracted from these 
fits are plotted as a function of temperature in Fig. [S] 

Both gaps are seen to decrease with temperature, 
falling to near zero at « 20 K. Above 20 K the fitted 
gaps are relatively constant and close to zero. 

III. DIFFRACTION MEASUREMENTS 
A. Experimental Details 

Neutron diffraction measurements were performed on 
the hot-neutron diffractometer D3 at the Institut Laue- 
Langevin. The instrument was used in unpolarized- 
neutron mode with a neutron wavelength of 0.84 A. The 
single crystal of Nao.ssCoO^ used for these measurements 
was cleaved from a rod grown in Oxford by the floating- 
zone method^ The crystal had a mass of 0.3 g and a 
mosaic spread of ~ 2 degrees. 

Magnetization measurements made using a SQUID 
magnetometer in Oxford confirmed the existence of the 
metamagnetic transition at ~ 9 T reported by Luo et al. 
at this composition^ 




5 10 15 20 25 

T(K) 

FIG. 5: (color online) The magnitudes of the two gaps as a 
function of temperature. Data points calculated from Eqn[I] 
using the values of J c , D and E derived from fits such as those 
in Fig|3 Solid curves are guides for the eye. Error bars were 
estimated by varying the two gap energies separately until the 
fit was no longer acceptable. 



The crystal was pre-aligned on the neutron Laue 
diffractometer Orient Express at the ILL and mounted 
on an aluminium pin using ceramic glue. The ideal setup 
for this experiment would be to align the c axis vertically, 
applying a vertical field, with the incident and scattered 
beams inclined at the Bragg angle to the horizontal in or- 
der to access the (0, 0, 3) reflection. An alternative setup, 
with a fixed horizontal incident beam, would require the 
cryomagnet holding the sample to be tilted by the Bragg 
angle (9b), with the detector lifted out of the plane (by 
29b), allowing access to the (0,0,3) reflection while still 
applying the field directly along the c axis. 

However, on D3 we were restricted to using a horizon- 
tal incident beam (ruling out the first setup), and also 
unable to tilt the cryomagnet (ruling out the second) . In 
order to access the (0, 0, 3) reflection we tilted the crys- 
tal c axis 7 degrees away from vertical, corresponding to 
the (0,0,3) Bragg angle for 0.84 A neutrons, and lifted 
the detector out of the horizontal plane by 14 deg. The 
10 Tesla vertical-field cryomagnet in which the crystal 
was mounted then allowed application of a field almost 
parallel to the c axis (though actually 7 degrees away 
from it). The field-induced (0,0,1) magnetic Bragg peak 
is expected to be larger than that at (0,0,3) due to the 
magnetic form factor, but with 0.84 A neutrons the scat- 
tering angle for (0,0,1) is too small to access with our 
setup. 



B. Results 

Figure shows the main results of the diffraction stud- 
ies of Nao.85Co02- All the measurements were made at 
Q = (0,0,3), scanning either field or temperature with 
the other external variable fixed. For ideal Naa:Co02 no 
structural Bragg peak is allowed at this position and no 
magnetic Bragg peak is allowed if the moments point 
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FIG. 6: Diffraction studies of Nao.ssCoO^. (a) Field scans 
at Q = (0,0,3) at constant temperatures of T = 1.8, 10 
and 30 K. Data at T = 10 K and 30 K are shifted up by 
500 and 1000 counts respectively, (b) Temperature scans at 
Q = (0, 0, 3), with zero applied field (open circles), and with 
H =9.6 T applied at 7 degrees to the c-axis (filled circles). 



FIG. 7: (color online) (a) The ordered A-type antiferromag- 
netic structure in the a-c plane, with an external magnetic 
field B applied parallel to the magnetic moments, (b) Above a 
critical field B s f the system undergoes a spin-flop transition, 
to a phase with spins at an angle 6 = (j> to the c-axis. 



(1) the field at which the spin-flop transition occurs shifts 
up slowly with and (2) the ordered magnetic mo- 
ment /X decreases with T. Both effects would result in 
the reduction of intensity with increasing temperature. 
However, the effect of (1) is too small to explain the dis- 
appearance of the signal by 30 K (see Luo et alM>), so 
we deduce that the signal disappears due to the reduc- 
tion of the magnetic moment to zero above the magnetic 
ordering temperature T m . 

In Fig. 03 we plot the temperature dependence of the 
scattering at Q = (0, 0, 3), in both zero applied field and 
9.6 T. We confirm that the signal induced by application 
of the magnetic field decreases to zero as the temperature 
is raised to ~ 20 K. 



C. Analysis 



along the c axis. 

Figure|SJi shows three field scans at Q = (0, 0, 3), per- 
formed at constant temperatures of T = 1.8, 10 and 30 K. 
The data at 10 K and 1.8 K have been shifted up by 500 
and 1000 counts respectively for clarity. 

At T — 1.8 K there is a large increase in the intensity 
of scattering at (0, 0, 3) between ^6 T and 9 T as the 
field increases, with the intensity appearing to flatten 
off between 9 and 10 T. At T = 10 K the increase in 
intensity has shifted up in field to start at ^8 T, and 
at T = 30 K the intensity remains constant with field. 
The increase in intensity with field at 1.8 K is consistent 
with the expected spin-flop transition because once the 
spins have rotated away from the c direction magnetic 
Bragg scattering is allowed at (0, 0, 3) provided that the 
ordering wavevector remains (0,0, 1). 

With increasing temperature two effects are at work: 



Figure |7| shows the magnetic structure of the low-field 
A-type antiferromagnetic (AF) phase (a), along with that 
of the spin-flop (SF) phase (b). As discussed above, the 
large increase in Bragg intensity at Q = (0, 0, 3) that 
occurs between ~ 6 T and 9 T (Fig. Efa)) appears to 
support the idea that a phase transition from an AF to 
a SF phase occurs at this field. To connect the spin-flop 
transition with the anisotropy gap described in Section 
II we extend the spin-wave model to include an external 
magnetic field applied along the c axis. We begin by 
calculating the effect of the field on the spin- wave modes 
of the AF phase. This enables us to calculate the critical 
field at which a phase transition would be expected to 
occur. A term is added to the original Hamiltonian to 
represent a vertical applied magnetic field, giving a new 
Hamiltonian: 

H' = H + g[i B Bj2s? , (?) 
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where H is the original Hamiltonian (Eqn|TJ, B is the 
magnitude of the applied magnetic field, and we assume 
.9-2. 

The spin-wave dispersion was derived from W as be- 
fore. The field further splits the two modes, i.e. the lower 
mode moves lower in energy, and the higher mode moves 
higher in energy, as the field is increased. The magnitude 
of the gaps at the magnetic zone center as a function of 
field, B, are then given by 



hwf a JB) = 2S p 2 + {2J C D + D 2 - E 2 



±2 V / /3 2 D(2J C + D) + J*E 2 



(8) 



where (3 — gnsB/2S. 

As the field is increased to a critical field, B c \ 1 a local 
instability occurs when the energy of the lower mode falls 
to zero and then becomes imaginary, and the system can 
no longer remain in the low-field AF phase (Fig. EH)- The 
critical field is determined by setting Ttojg ap (B c {) = 0, 
resulting in the expression 



B c i — 



2S_ 



^(D-E)(2J C + D-E) 



(9) 



B c i is the critical field at which we would expect the 
system to 'flop' out of the A-type antiferromagnetic phase 
(Fig. Uk) into the spin-flop phase (Fig. 0d) based on the 
closing of the gap with increasing field. 

A similar calculation can be performed by considering 
the spin- wave modes in the SF phase, shown in Fig[7fb), 
with 9 = <f>. In this phase there are two modes at high 
field, and as the field is decreased there exists a critical 
field B C 2 at which the lower mode vanishes and the SF 
phase is no longer stable. At this field, B C 2, the system 
returns to the AF phase. By calculating the spin-wave 
dispersion in the SF phase we find the angle at which the 
spins lie in the a-c plane, as a function of field, 



9 = arccos 



[iB 



2S 2 (2J C -D + E)J ' 
and an expression for the critical field, B C 2- 



(10) 



b c2 =^J {d -3 2 ±-°z e)2 . (id 



9l L B 



(2 J c + D-3E) 



In a spin- flop transition B c \ > B C 2, so hysterisis is 
possible. By evaluating EansPandlTTl using the values 
for the exchange and anisotropy constants determined by 
fitting the inelastic data on Nao.7sCo02 at 1.5 K (Sec- 
tion II C), we find values of B cl = 8.19 ± 0.2 T and 
B c i = 8.15±0.2 T. Therefore the hysterisis is expected to 
be too small to measure, and we take B c = 8.2 ±0.2 T as 
the critical field. This value is in strong agreement with 
the experimental data shown in Fig. 01. When the spins 
'flop' into the SF phase we calculate 9 = 87.8±0.5 degrees 
(from Eqn. HOfl , so the spins lie almost aniferromagnet- 
ically perpendicular to the applied magnetic field. We 



note that a mean-field approximation to this calculation 
also gives values B c = 8.2 ± 0.2 T and 9 « 88 degrees, 
and that corrections to this calculation to account for the 
field offset of 7 degrees from vertical are small and within 
the given errors. 

An estimate of the magnetic moment was derived from 
the diffraction data by comparing the integrated inten- 
sities of a set of nuclear reflections with the integrated 
intensity of the (0,0,3) magnetic Bragg peak, using the 
A-type antiferromagnetic model for the magnetic struc- 
ture. From this we estimate a magnetic moment of 
0.11±0.07/ib, m agreement with values derived by other 
techniques M> For this calculation we used the free-ion 
form factor for Co 4+ , which may not be a good approxi- 
mation given the metallic nature of Nao.85Co02- 



IV. DISCUSSION 

We have seen that the spin-wave dispersion does con- 
tain a clear gap, and that at low temperature (1.5 K) it is 
well described by our simple Hamiltonian extended to in- 
clude two anisotropy terms (Eqn. 2J. While including a 
term to describe a uniaxial anisotropy (D) seems logical 
given that the spins do lie along the c axis, the need to 
introduce also an in-plane anisotropy term, with parame- 
ter E nearly two thirds as large as the uniaxial parameter 
D, was unexpected. Although the form of the in-plane 
anisotropy term may need more careful thought, we have 
shown that in-plane anisotropy clearly exists in the com- 
pound: the data cannot be described with the uniaxial 
anisotropy alone. The magnitudes of the anisotropics D 
and E (0.096 and 0.059 meV respectively at 1.5 K) are 
very small in comparison to the exchange parameters J c 
and J a b (12.2 and 6 meV). To refine the model further 
with a more realistic form for the spin anisotropy would 
require a more detailed consideration of the source of the 
anisotropy. 

By fitting the same model to data measured at various 
temperatures up to 24.4 K we have shown that both gaps 
disappear quite suddenly at approximately 20 K. From 
the proximity of this temperature to the bulk magnetic 
transition at T m « 22 K found in magnetization studies 
we infer that the spin waves have the same origin as the 
T m transition. This may seem obvious, but the strength 
of the spin-wave scattering in comparison to the small 
size of the ordered moment made it important to confirm 
that the observed spin excitations were really associated 
with the magnetic order. 

We have successfully observed the magnetic Bragg 
peak at Q = (0, 0, 3) in Nao.85Co02 with unpolarized 
neutrons by inducing a spin-flop transition in a verti- 
cal magnetic field. The transition occurred at a field of 
B s f « 8T, but was broad, with a width of approximately 
3 T. It is possible that the broadening of the transition 
may be due to disorder in the structure, which would lead 
to a spread of values of the interlayer exchange constant 
J c and in turn lead to a range of values for B s f . 
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We have presented a simple calculation of the spin- 
flop transition field that would be expected based on the 
anisotropy parameters derived from the measurement of 
the spin gap in Nao.7sCo02- The calculated value for 
this critical field of 8.2 ± 0.2 T is in strong agreement 
with the observed value for Nao.ssCoCb of B s f w 8 T, 
We take this as further evidence that the A-type anti- 
ferromagnetic structure with its associated gapped spin 
excitations is the same phase marked by the transition 
temperature T m «22 K in magnetization data. 

Although we have confirmed the magnetic structure 
and characterized the excitation gap, the nature of the 
magnetic ground state in Na^OoC^ is still in question. 
In an ionic picture there are (1 — x) Co 4+ ions carry- 
ing spin S = \ in a background of x Co 3+ ions, which 
are usually assumed to be non-magnetic. However, it is 
uncertain whether these ions would be ordered or clus- 
tered or randomly distributed, 8 and there is also evidence 
that Na^CoO^ is a good metal, which would suggest 
that an itinerant picture might be more appropriated 
A weakly itinerant ground state with strong spin fluctu- 
ations would be consistent with both the small ordered 
moment (< 0.2/i#) and the fact that the energy scale of 
the magnetic excitations is much greater than the mag- 
netic ordering temperature. 

In this vein, the clean gap we have observed in the 
spin- wave dispersion points towards a system of local mo- 
ments with a small symmetry-breaking anisotropy field, 
which is supported by our successful description of the 
data using a simple spin-wave model based on localized 
Co spins. However, the large intrinsic widths are not con- 
sistent with a model of purely localized Co ions, and can 
be taken as evidence of more metallic behavior. These 
features of the magnetism of Na x Co02 need to be taken 
into consideration when assessing theories of the metallic 
state of this material. 



V. CONCLUSIONS 

We have studied the low-energy region of the excita- 
tions in Nao.75Co02 and revealed that the excitations are 
gapped. At low temperature the data were modelled best 
with a simple spin-wave model with two non-degenerate 
modes, providing evidence for both uniaxial and in-plane 
anisotropy in the compound. Gaps from both modes were 
found to close at the bulk magnetization temperature 
T m w 22 K. We therefore conclude that the spin waves 
observed in Nao.7sCo02 are associated with the magnetic 
phase transition seen previously in magnetization data. 

Furthermore, we have confirmed the A-type an- 
tiferromagnetic structure underlying these excitations 
by performing neutron diffraction measurements on 
Nao.85Co02 in a large vertical applied magnetic field. A 
spin-flop transition occurs at « 8 T allowing measure- 
ment of the magnetic Bragg peak at Q = (0,0,3). The 
results support the assumption that the spins lie along 
the c axis. We have calculated the field needed to over- 
come the anisotropy associated with the gap in the spin- 
wave dispersion and shown that it is in good agreement 
with the observed spin-flop transition field. 
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